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National Academy of Science 
P. E. Schmid 

NASA /Goddard Space Flight Center 
ABSTRACT 

This note deals with the measurement of the total electron content 
of the ionosphere at the Goddard Space Flight Center, looking towards 
the geostationary satellite ATS-3 during the solar eclipse of March 7, 
1970. Obscuration at this site was near total. Faraday rotation was 
measured using a stationary circularly polarized antenna and a dual 
channel phase lock receiver tuned to 137, 350 MHz. By comparing the 
electrical phase of the two opposite circularly polarized components 
a continuous chart recording was made of Faraday rotation versus local 
time. 

A depletion of about 25% in electron content was observed from 
first contact to the time of minimum electron content. The time 
variations of the electron content during the eclipse are briefly ex- 
amined in the light of current theories of ionospheric processes. 


CONTENTS 


Page 

1. INTRODUCTION - THE MARCH 7, 1970, SOLAR ECLIPSE. ... 1 

2. FARADAY ROTATION MEASUREMENT . 1 

3. COMPUTATION OF TOTAL ELECTRON CONTENT 2 

4. VARIATIONS OF THE TOTAL ELECTRON CONTENT (TEC) 

DURING THE ECLIPSE 5 

5. DISCUSSION OF THE ECLIPSE EFFECTS 6 

ACKNOWLEDGEMENT 9 

REFERENCES 9 

ILLUSTRATIONS 

Figure Page 

1 Solar Eclipse 10 

2 Polarization Measurement Scheme 11 

3 Faraday Rotation from ATS-3 . 12 

4 M -Factor Variation with Height along Ray Path Linking GSFC 

with ATS-3 13 

5 Total Electron Content and Optical Obscuration Function 

During the Solar Eclipse of March 7 , 1970 ......... 14 

2 

6 Relative Faraday Rotation Angle Versus (f 0 f 2 ) . 15 

7 Variation of SHMAX and SCAT During Eclipse (Raw Data) ... 16 

8 Smoothed Variatio n of SHMAX and SCAT along with TEC 

and Obscuration Function 17 

9 Variation of SCAT, RATIO, and HMAX During the 

Eclipse 18 

iv 



1. INTRODUCTION - THE MARCH 7, 1970, SOLAR ECLIPSE 

The lunar shadow associated with the March 7, 1970, solar eclipse passed 
along the east coast of the United States between 1315 and 1400 hours EST. The 
shadow path is indicated in Figure 1. During the eclipse the NASA geostationary 
spacecraft Applications Technology Satellite-3 (ATS-3) was positioned at 275° 

East longitude, 0.53° South latitude, and 42,130km altitude. Thus, the ATS-3 
137.350 MHz transmission as recorded at the Goddard Space Flight Center (GSFC) , 
Greenbelt, Maryland, traversed an ionosphere which during the eclipse experi- 
enced near total obscuration. The basic measurement made during the eclipse 
was relative spatial polarization of the linearly polarized ATS-3 transmissions 
as a function of time. 

2. FARADAY ROTATION MEASUREMENT 

As is well known, changes in the total electron content through the Earth’s 
ionosphere can be related directly to Faraday rotation measurements. At GSFC 
Faraday rotation of the ATS-3 VHF signals were recorded during the eclipse by 
means of a two channel phase lock receiver, each channel having as an input one 
of the two circular components derived from the incoming linear polarization. 
Figure 2 presents the basic measurement technique where the linearly polarized 
signal is resolved into right and left hand circular components and then multi- 
plied in an electronic phase detector whose full scale output is 2 7r radians where 
each n of rotation in spatial angle 0 results in one cycle of electrical phase out- 
put. The receiver output, as shown in Figure 3, was recorded directly by an 
Esterline Angus chart recorder running at a rate of 1. 5 inches per hour. 
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As shown in the records presented in Figure 3, during a normal day (such 
as February 25, 1970) ionization increases steadily during the morning hours up 
to the afternoon peak and then decays steadily as solar illumination decreases. 
Again with reference to Figure 3, a pronounced phase reversal occured during 
the eclipse as the lunar shadow intercepted the line-of-sight between the GSFC 
site and ATS-3 and ionization abruptly decreased. At approximately 1400 hours 
local time the integrated electron content again increased up to the usual diurnal 
reversal. 

3. COMPUTATION OF TOTAL ELECTRON CONTENT 

The Faraday rotation of the plane of polarization of a linearly polarized 
wave is given by the well known expression 

hs 

0.0297 r (1) 

Q = I N(h) M dh 

f 2 J 

o 

where M = H cos 6 sec i 

f = frequency of the wave 

N(h) = electron density at a height h in the ionosphere 
H = the geomagnetic field intensity in MKS units 

0 = the angle between the ray and the geomagnetic field at 

the height h 

1 = the angle between the ray and the vertical at the height h 
h s = vertical height of the satellite 
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The quantity M was computed by knowing the geometry of the ray path con- 
necting the receiver with the satellite and using the spherical harmonic expan- 
sion of Hendricks and Cain (1966) for the geomagnetic field. Figure 4 shows the 
variation of the quantity M with height along the ray path. It may be seen that 
M is a slowly varying function of height when compared with the possible varia- 
tion of the electron density with height in the ionosphere. Hence, it is possible 
to approximate equation (1) by 

hs 

0. 0297 _ ? (2) 

ft M j N(h) dh 

f 2 J 

o 

where M is the value of M at some ionospheric height chosen so as to make 
equations (1) and (2) numerically most equivalent. In the present calculations 
the value of M was chosen to be that of M at a height of 350 km along the ray 
path. 

In earlier studies of ionospheric electron content during a solar eclipse 
signals from near-Earth satellites had to be used and this led to the variation 
of M over the duration of the eclipse due to the motion of the satellite and the 
consequent variation of the ray path. This is true of eclipse studies by the lunar 
radar technique too, though to a lesser extent. In the case of the March 7, 1970, 
eclipse, however, it was possible to make continuous measurements of Faraday 
rotation over a constant ray path by using the geostationary satellite ATS-3. 
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The accumulated Faraday rotation of the plane of polarization was plotted 
against time by using an arbitrary reference angle at some time t(0844 EST) and 
adding observed rotations relative to it from the original recording shown in 
Figure 3. 

Figure 5 shows such a plot. For comparison the percentage obscuration of 
the solar disk as viewed from a point at a vertical height of 300 km on the ray 
path linking the GSFC receiving site with the ATS-3 is al£>o shown in the Figure. 

It is possible to compute changes in the total electron content (TEC) from 
the relative rotation plot by using equation (2). However, the computation of the 
absolute values of the TEC requires a knowledge of the absolute value of Faraday 
rotation of the plane of polarization. An attempt was made to resolve this ambi- 
guity in the absolute angle of rotation by comparing it with ground based meas- 
urements. It is well known that the TEC is proportional to the maximum electron 
density of the ionospheric F2 layer and this fact has been used by Nakata (1966) 
to deduce the value of the absolute angle of Faraday rotation. Therefore, the 
relative angle of rotation was plotted against the square of the critical frequency 
of the F2 layer as measured at Wallops Island, Virginia. This was done be- 
cause the maximum electron density, N m # is proportional to the square of the 
critical frequency of the F2 layer. Figure 6 shows the resulting plot. It may be 
seen that the plot may be approximated by linear segments of differing gradients. 
The segment nearest to the vertical axis was fitted with a best fitting straight 
line, using the least squares criterion. This line is also shown in Figure 6. 
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The intercept of this straight line with the Y axis was taken to be the reference 
angle of rotation, The segment nearest to the vertical axis was chosen because 
it is numerically sound to assume a linear gradient only over a short range of 
values and the segment chosen needs the least amount of extrapolation to inter- 
cept the Y-axis. 

Once the reference angle of rotation is known, it is easy to compute the 
absolute angles of rotation and hence to deduce the values of TEC using equation 
(2) . The results of such a calibration are shown on the right hand side of 
Figure 5. 

4. VARIATIONS OF THE TOTAL ELECTRON CONTENT (TEC) DURING 

THE ECLIPSE 

The TEC variations show the normal diurnal increase around noon at t s 
beginning of the eclipse. This increase continued for nearly 25 minutes after 
first contact. At 1245 EST the effect of the eclipse became apparent by a sharp 
decrease which continued until the minimum during the eclipse was reached at 
1356 EST. The recovery of the TEC during the second phase of the eclipse was 
at first quite slow giving the impression of a flat minimum for about 30 minutes. 
From about 1425 EST the TEC started to increase rapidly until about 1510 EST 
when the post -diurnal -maximum trend took over and the TEC started to decrease 
again. 

It may be observed from Figure 5 that the minimum in the TEC lags behind 
the time of maximum obscuration (1340 EST) by 16 minutes. The TEC decreased 
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from a value of about 6 x 10 17 el/m 2 to abort 4.4 x 10 17 el/m 2 during the eclipse. 
Thus , the depletion in the columnar electron content during the eclipse was about 
1,6 x 10 17 el/m 2 or 25%. 

5. DISCUSSION OF THE ECLIPSE EFFECTS 

It is interesting to compare the variations of TEC with other ionospheric 
parameters like the electron content below the level of maximum ionization 
density (SHMAX) and the thickness paramete* of the F2 layer (SCAT) as these 
yield more information about the layer structure. During the March 7, 1970, 
eclipse, intensive ground based sounding of the ionosphere was made at Wallops 
Island, Virginia, which was in the track of the eclipse. The analysis of these 
ionograms were done by the ESSA group to deduce parameters like the electron 
density at different heights N(h), SHMAX, SCAT, etc. Using the data furnished 
by ESSA, the time variations of SHMAX and SCAT were plotted and this is shown 
in Figure 7. It may be seen from the Figure that the bottom side electron con- 
tent was reduced by nearly 40% during the eclipse. The value of SCAT which can 
be considered as the scale height at the height of the maximum electron density, 
HMAX, remained near 50 km for most of the time except near 1350 EST when 
there were rather large fluctuations between about 72 km to 36 km. 

In fact, it may be seen from Figure 7 that the SHMAX and SCAT data show 
considerable oscillation and in the case of SCAT the oscillations make it difficult 
to see the trend of the main variations. An attempt was made to smooth the data 
by fitting least-squares polynomials to them. Figure 8 shows 7th degree poly- 
nomial fits to the SHMAX and SCAT data along with the TEC variations. It nmy 
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be observed from the Figure that the bottom side electron content more or less 
follows the total electron content. However, its response to the eclipse is faster 
in both phases of the eclipse. The smoothed variation of the SHMAX data shows 
a delay of about 8 minutes between the minimum and the time of maximum ob- 
scuration. The smoothed SCAT data show two minima with a maximum around 
the time of the minimum in TEC. 

The ratio of the TEC to SHMAX is another parameter defining the layer 
structure. This quantity which we shall call RATIO is plotted in Figure 9 along 
with SCAT and HMAX, the height of the maximum electron density of the F2 
layer. It may be seen that the value of RATIO during the eclipse fluctuates be- 
tween 2. 8 to 3. 6. The value of RATIO for the classical Chapman layer is 3. 15. 
It is interesting to see that while the SCAT variations are roughly in phase with 
those of HMAX, the RATIO variations are roughly anti-phase to those of SCAT. 
This is what one would normally expect in view of the fact that HMAX divides 
the top side from the bottom side and an increase of HMAX will normally de- 
crease RATIO. This is particularly true in view of the fact that SCAT is related 
to the electron and ion temperatures (Stubbe, 1970) and will normally increase 
with increase of HMAX. An increase of SCAT will further increase SHMAX, 
However, the effects of downward diffusion of ionization coupled with the cooling 
of the plasma during the eclipse may lead to departures from this simple pattern. 

During the first phase of the eclipse RATIO is seen to increase up to about 
1320 EST. This is probably due to the fact that the bottom side ionosphere is 


responding more rapidly than the top aide to the eclipse. The eventual reduc- 
tion of the top side ionization coupled with the downward diffusion of ionization 
(Stubbe, 1970) reduces RATIO later on. After the time of maximum obscuration 
the bottom side ionosphere recovers more rapidly than th(* top side which makes 
RATIO decrease still further. Finally the response of tho top side ionosphere 
to the second phase of the eclipse coupled with the decrease in downward diffu- 
sion of ionization (Stubbe, 1970) leads to the increase of RATIO again. 

Stubbe (1970) has made a recent theoretical study of the time variations of 
the ionosphere over Wallops Island during the eclipse of March 7, 1970. He did 
this by the numerical solution of ten coupled differential equations , including 
the continuity equations for different ions , the heat conduction equations and the 
equations of motion describing neutral winds. The time dependent solutions 
yield such quantities as electron and ion concentrations and temperatures and 
hence layer parameters such as integrated electron content and semithickness 
of the layer. Stubbe predicted the following for the March 7, 1970, eclipse: 

1. The ratio of the TEC to the bottom side electron content will have a 
sharp maximum of 3. 3 at the time of the maximum obscuration as 
compared to the value 2. 6 for the reference atmosphere. 

2. The TEC will have a minimum at the time of maximum obscuration 
which is about a factor of two smaller than the corresponding value 


in the normal ionosphere. 


\\ 


From our study we note that RATIO increased from a value of about 2. 9 at 
1200 EST to a maximum of 3. 6 at 1320 EST. These values are in good agree- 
ment with Stubbed predictions. The TEC has a minimum during the eclipse 
which is smaller by a factor of 1. 5 than the corresponding value in the normal 
ionosphere. However, as pointed out earlier, this minimum in TEC is not 
coincident with the time of maximum obscuration but lags behind it by 16 minutes. 
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FIGURE CAPTIONS 


Figure 1. Solar Eclipse 

Figure 2. Polarization Measurement Scheme 

Figure 3. Faraday Rotation from ATS-3 

Figure 4. M-Factor Variation with Height along Ray Path Linking GSFC with 
ATS-3 

Figure 5. Total Electron Content and Optical Obscuration Function During 
the Solar Eclipse of March 7, 1970 
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Figure 6. Relative Faraday Rotation Angle Versus (f 0 f 2 ) 

Figure 7. Variation of SHMAX and SCAT During Eclipse (Raw Data) 

Figure 8. Smoothed Variation of SHMAX and SCAT along with TEC and 
Obscuration Function 

Figure 9. Variation of SCAT, RATIO, and HMAX During the Eclipse 
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Figure 1 . Solar Eclipse 
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Figure 2. Polarization Measurement Scheme 
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FARADAY ROTATION FROM ATS-3 
SOLAR ECLIPSE OF 7 MARCH 1970 
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Figure 3. Faraday Rotation from ATS-3 
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Figure 4. M-Factor Variation with Height along Ray Path 
Linking GSFC with ATS-3 
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Figure 5. Total Electron Consent and Optical Obscuration Function During 
the Solar Eclipse of March 7, 1970 
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Figure 9. Variation of SCAT, RATIO, and HMAX during the Eclipse 
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